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Abstract 

A method for homogenization of a heterogeneous (finite or periodic) elastic compos- 
ite is presented. It allows direct, consistent, and accurate evaluation of the averaged 
overall frequency-dependent dynamic material constitutive relations. It is shown 
that when the spatial variation of the field variables is restricted by a Bloch-form 
(Floquet-form) periodicity, then these relations together with the overall conserva- 
tion and kinematical equations accurately yield the displacement or stress mode- 
shapes and, necessarily, the dispersion relations. It also gives as a matter of course 
point-wise solution of the elasto-dynamic field equations, to any desired degree of 
accuracy. The resulting overall dynamic constitutive relations however, are general 
and need not be restricted by the Bloch-form periodicity. The formulation is based 
on micro-mechanical modeling of a representative unit cell of the composite pro- 
posed by Nemat-Nasser and coworkers; see, e.g., [1] and [2]. 

We show that, for a micro-structured elastic composite, the overall effective mass- 
density and compliance (stiffness) are always real-valued and positive, whether or 
not the corresponding unit cell (representative volume element used as a unit cell) 
is geometrically and/or materially symmetric. The average strain and linear mo- 
mentum are however couple and the coupling constitutive parameters are always 
each others complex conjugates for any heterogeneous elastic unit cell, such that 
the overall energy-density is always real and positive. In this paper, we have sought 
to separate the overall constitutive relations which should depend only on the com- 
position and structure of the unit cell, from the overall field equations which should 
hold for any elastic composite; i.e., we use only the local field equations and material 
properties to deduce the overall constitutive relations. 

It is shown, by way of an example of a bi-layered composite, that dispersion 
curves obtained by our method accurately produce the exact results of Rytov [3]. 
The method is also used to calculate the effective parameters for a 2-layered com- 
posite and the results are compared with those of homogenization based on the 
field integration of the exact solution (Willis [4], and Nemat-Nasser et al. [5]), and 
certain relevant issues are clarified. Finally the method is used to homogenize both 
a symmetric and a non-symmetric 4-layered composite and the results for the sym- 
metric case are compared with those reported by Nemat-Nasser et al. [5] as well 
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as the exact solution. Thus, this method provides a powerful solution and homoge- 
nization tool to use in many cases where the unit cell contains inclusions of complex 
geometry. 

Key words: Homogenization, Microstructured Elastic Composites 
PACS: *43.20.Gp, *43.20.Jr, 62.20.D- 



1 Introduction 



The recent interest in the characterization of the overall dynamic properties 
of composites with tailored microstructure, necessitates a systematic homog- 
enization procedure to express the dynamic response of an elastic composite 
in terms of its average effective compliance and density. The resulting homog- 
enized parameters, once restricted by a Bloch-form periodicity, must give the 
composite's dispersion relations as well as the corresponding strain or stress 
mode-shapes. 

In the present paper, we propose a method for homogenization of microstruc- 
tured elastic composites that directly provides the overall frequency-dependent 
dynamic material parameters. In addition, when the spatial variation of the 
overall field quantities is restricted to follow a Bloch-type periodicity, then the 
overall conservation and kinematical relations yield the corresponding stress 
and displacement mode-shapes and dispersion relations. The method does not 
require the values of the field variables within the unit cell. Indeed, it actually 
gives as a matter of course the point-wise solution of the elasto-dynamic field 
equations to any desired degree of accuracy. The method is inspired by the 
micro-mechanical homogenization of the static properties of a typical unit cell 
of a composite, originally proposed by Nemat-Nasser and Taya [6], and fur- 
ther developed in [7] and [2] . A similar homogenization method has been used 
by Amirkhizi and Nemat-Nasser [8] to calculate the effective electromagnetic 
properties of a periodic composite. Here however, we clearly distinguish the 
overall constitutive relations which should reflect the composition of the unit 
cell, from the overall conservation and kinematical relations which should hold 
for any (homogenized) elastic composite with any unit cell. In other words, 
while we do use the local (at each point within the unit cell) material proper- 
ties and do enforce the conservation laws at the local level to obtain the overall 
constitutive relations, the overall conservation and kinematical relations are 
not used to derive these overall constitutive relations. Therefore, the overall 
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constitutive relations can be used to solve initial-boundary value problems 
for waves of wavelengths that are large relative to the micro-structural length 
scale. 

Homogenization for calculating the overall properties of composites, based 
on the integration of the field variables, has been proposed by a number of 
researchers. For electromagnetic waves, see, for example, Pendry and Smith 
[QlllOfll] . Amirkhizi and Nemat-Nasser [12], Bensoussan [13], and Sihvola [14] . 
For elasto-dynamic waves, Willis [1] has presented a homogenization method 
based on an ensemble averaging technique of the 'Bloch' reduced form of the 
wave propagating in a periodic composite. An analogous approach for the 
calculation of effective compliance and mass density of a layered composite 
was worked out in some details by Nemat-Nasser et. al. p] who showed that 
the corresponding results yield the exact dispersion relations. 

In what follows, we outline our approach in terms of a layered composite. We 
compare, by way of an example of a 2-layered composite, the corresponding 
dispersion results with those obtained using the exact solution. Their close cor- 
respondence indicates that our micro-structurally-based method may be used 
to calculate the dispersion relation for more complex 2- and 3-dimensional 
cases, where exact solutions are not available, as we shall report elsewhere. 
Here we also use our approach to homogenize a 4-layered composite and com- 
pare the results with those obtained by the field integration of the stress and 
displacement mode-shapes obtained using the mixed variational formulation of 
Nemat-Nasser et al. [T5|16fl7] . It is shown that the micro-structurally-based 
method gives homogenized results which converge to the field integration- 
based homogenization results that are based on either the exact solution or 
the mixed variational formulation. Thus, the micro-structurally-based method 
can be used to evaluate the effective dynamic parameters of complex unit cells, 
without a priori requiring explicit point-wise values of the field variables, but 
instead producing these quantities if desired. 



2 Micro-structural homogenization of layered composites 

Here we present a homogenization method based on a micro-mechanical con- 
sideration of the volume averages of the field variables, viewed as measurable 
macroscopic physical quantities. We express the solution to the elasto-dynamic 
equations of motion as the sum of the volume average and a perturbation due 
to the heterogeneous composition of the unite cell, 

q = qo + qp (1) 
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where Q represents any of the field variables, stress, a, or velocity, u. The aim 
is to derive a set of constitutive relations for the overall averaged parts of the 
field variables, using the local elasto-dynamic equations of motion and con- 
stitutive relations. This then provides the homogenized frequency-dependent 
material parameters. In what follows, we describe our approach using a peri- 
odic composite, but the final constitutive relations also apply to a finite unit 
cell. 

Consider harmonic waves in an unbounded elastic composite consisting of a 
collection of bonded, identical unit cells, = {x : —a/2 < x < a/2}, which 
are repeated in the x-direction, and hence constitute a periodic structure. In 
view of the periodicity of the composite, we have p{x) = p{x + m a) and 
C{x) = C{x + m a); here m is an integer, p{x) is the density and C{x) is 
the modulus of elasticity. For time harmonic waves with frequency u, the field 
quantities are proportional to e^*"^*. For waves with arbitrary wavenumber q 
and unrelated frequency u, we express the field variables in the following form: 



F{x, t) = Re [F{x) exp[i{qx - ut)]] (2) 

where F represents the field variables, stress, a, strain, e, momentum, p, dis- 
placement, u, or velocity, m, whereas F represents their periodic parts (a, e, p, 
u, u). The representation, Eq. separates the time harmonic and macroscopic 
factor from the microscopic part of the field variables. Even for a finite unit 
cell, the Fourier series solution of the microscopic part is periodic so that this 
solution satisfies, F{x) = F{x + m a). We emphasize here that the frequency 
and the wavenumber, u and q, are, at this point, unrelated and arbitrary. 

The local conservation and kinematic relations are, 

Vcr = —iup 

(3) 

Vn = —ico£ 



where V — )■ ^ + ig. The corresponding constitutive relations are. 



e = D(x)a 

^ ' (4) 
p = p{x)u 

where D{x) is the compliance, the inverse of the elastic modulus, and p{x) is 
the mass-density. These local material parameters represent the structure and 
composition of the unit cell. 
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Now we replace the heterogeneous unit cell with a homogeneous one hav- 
ing uniform density po and compliance Dq. In order to reproduce the strain 
and momentum of the actual unit cell, field variables eigenstress, and 

eigenvelocity, U{x), are introduced. These quantities are then calculated, us- 
ing the basic local field equations and constitutive relations. The idea stems 
from the polarization stress or strain that was originally proposed by Hashin 
[T9] and further developed by Hashin and Shtriktman [20||2T] and later by 
others, in order to construct energy-based bounds for the composite's overall 
elastic moduli. The basic tool in these works has been the result obtained by 
Eshelby [22j in three dimensions and earlier by Hardiman [23| in two dimen- 
sions, that the stress and strain are constant within an ellipsoidal (elliptical in 
two dimensions) region of an infinitely extended uniform elastic medium when 
that region undergoes a uniform transformation corresponding to a uniform 
inelastic strain. 

Here, however, we present a different tool that can be used to actually calculate 
the point-wise values of the elasto-dynamic field variables to any desired degree 
of accuracy. For this, we require that the actual values of the field variables at 
every point within the homogenized and the original heterogeneous unit cell 
be exactly the same. To ensure this, we require that the following consistency 
conditions hold at every point within the unit cell: 

e = Da = Do{a~t) 
p = pu = po{u - U) 

The eigenstress and eigenvelocity fields are zero in regions where the mate- 
rial properties of the heterogeneous unit cell are equal to the chosen uniform 
material properties, Dq and po- From Eqs. ([3l [5]) we have, 

VV + z/V = z/^E - ^^V^ (6) 
V^u + v'^ii = u^U VS (7) 

lUjpo 

where z/^ = u^pqDq. We now consider a Fourier series solution of the above 
equations. With ^ = ±2mi/a,n ^ 0, we set, 

F{x) = F' + FP= {F)n + E HO^'^" (8) 
5 



(5) 



{F)n = ^J^F{x)dx 



(9) 



F(n = 1-1 F(x)e-'^''dx 



(10) 



where {F)q represents the averaged value of the field variable, F, over the unit 
cell, and F{C,) represents the Fourier coefficient of the corresponding local per- 
turbation due to heterogeneity. From Eqs. (El [7j) we now obtain the following 
Fourier coefficients for the stress and velocity fields: 



^(0 



- (^ + 



■m - 



UjD,[p^-{i + qy\ 



'111 



u{0 



- (^ + qy- 



[12) 



Therefore, the stress and velocity fields can be expressed as a sum of their 
average and their periodic components. 



a{x) = (a) + 2^ e 



loDq Jn 



(13) 



U[X) 



^(0^ / (^iy)e-'^'dy-^^ I t{y)e-^^ydy 
\l Jn ujpo ii Jn 



(14) 



^(0 



[u^ - (e + q?] 



K - (e + qf] 



(15) 



where (tt) and (cr) are the average values of the velocity and stress fields, 
respectively, taken over a unit cell. Note that -8(0 has the form. 



B{0 = iA{0 + qA{0 



(16) 



To make the homogenized unit cell point- wise equivalent to the original hetero- 
geneous unit cell, the homogenizing fields are required to satisfy the following 
consistency conditions: 
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D{x)[{cr)+aP] = Do[{a) + aP-t] 



(17) 



p{x)[{u) + uP]=Po[{u)+uP-U] (18) 

The periodic parts of the velocity and stress fields, from Eqs. (fl^ [T3|) . are 
now substituted into the above equations. This gives a set of 2 coupled in- 
tegral equations which yields the required homogenizing stress and velocity 
fields that exactly and fully replace the heterogeneity in the original medium. 
Our immediate concern, however, is not the point-wise representation of the 
heterogeneous medium, but, rather, it is the determination of the averaged 
field values, though the solution technique also yields the point-wise values of 
the field values as well. 

We now divide the unit cell into a number of subregions, fla- Then, we average 
the perturbation fields over each such subregion to obtain. 



1 



a^{x)dx 



(19) 



= E9a{0 {MO^ I (^iy)e-'^'dy-^^ f t{y)e-^^ydy] 

9a{0 = 7^ / ^'^""dx (21) 

(y. 

We now replace the integrals in Eqs. (IT^ [2U]) by their equivalent finite sums 
and set, 



^ jj{y)e-^^ydy ^Y.fm{-i)F, 



^ (22) 



n 

Fp = {F)n, 

Eqs. (fT9l [20ll then yield the following expressions: 
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a 







(23) 



ujpo 



(24) 



where the repeated index, [5, is summed over the number of subregions, [5 
1, . . . , a. The coefficient matrices in the above equations are defined by, 



(25) 



In these equations, /3 is not summed. Averaging the consistency conditions 
over each subregion a and using Eqs. ( I23l 12^ . we have. 



-BapUp 



(26) 



-Bap^l3 



^al3 H 

Pa - PO 



(27) 



Aal3 — faAaj3 

Bap = faBap o. uot summed 



(2J 



Note that Aq,^ and Baj3 are geometric quantities, independent of the material 
properties of the unit celL Expressions (!26| |27|) are 2a hnear equations which 

can be solved for a number of and a number of Up in terms of the average 
stess (a) and average velocity {u). We express the solution in the following 
matrix form: 



{S} = {$}(a) + — {*}(^) 
. 1 

{u} = -{e}(a) + {r}(^) 

Po 



(29) 



where, 
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{*} 



[Ao] + ^[B] [A,]-^ [B] 



[Ad] + -.[B] [A,]-^ [B] 



{/} 



[A,] + - [B] [Ad 



-1 



Bl 



-1 



[B] [A,]-'{f} 



\B] \A 



1-1 



D 



{/} 



{r} = [- [A, 

[AJ . = Aap + 



[Bl \Au]-' [Bl 



'■PiafS 
I T 



-r faDo 

Oal3 

Pa - PO 



{/} 



{fV ={/l,/2,...A} 



(30) 



Now we average the consistency conditions over the entire unit cell to express 
the average strain and average momentum in terms of the average stress and 
average velocity. Noting that the average of the periodic parts vanish when 



taken over the entire unit cell, we have, 



{e) = Do[{a) - (S)] 
= D{a) + Si{u) 



(31) 



ip) = Pom - m 

= S2{a) + p{u) 



(32) 



where 



D = Do[l - {ffW] 



Si 

S2 



Do 

-poifV-m = 

Po 



p = po[i-{fym 



(33) 



Eqs. fl^ 132]) are our final constitutive relations for the homogenized com- 
posite. In these equations, D and p are the overall homogenized composite's 
compliance and mass-density, respectively. As is shown below, D and p are 
alway real-valued and positive. The coupling terms, Si and §2, are always 
each others' complex conjugate. These attributes hold whether the unit cell is 
symmetric or non- symmetric. The coupling terms of course vanish if the unit 
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cell is homogeneous. They are real-valued and equal when the unit cell has 
reflective symmetry. 



2. 1 Mathematical Structure of Constitutive Relations 



We note that the matrices [Ad], [Ap], and [B], given by, 



[Ad],^ = E [U9.{i)fp9p{-i)A{i) + f^g^{-S)f,gp{i)A{-i)] 



5>0 Pa - PO 



a/3 



(34) 



are Hermitian for real values of the wavenumber. For any Hermitian matrix 
[M], one has. 



[[M]"']* = [[M]*]"' = [M]"' 



where asterix denotes a Hermitian transpose. Using this and the basic proper- 
ties of transpose of products of matrices, it follows that both — [Ad] + ^ [B] [Ap] 



and 



- [Ap] + ^ [B] [Ad]-^ [B 



-1 



are Hermitian matrices. Moreover, for a 
Hermitian matrix [M], it can be shown that the scalar {/}^ [M] {/} is always 
real- valued. Therefore, from Eq. f l30p . 



D = D^ 
P = Po 



1 - {fV 
I -iff 



-[Ad] + 47[B] [Ap]-^ [B 



AJ + ^ [B] [Ad]-' [B 



Pi ' 1 



-1 



{/} 
{/} 



(35) 



are always real-valued. We also have the following identities: 
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1 

UJ 

1 

u 

1 

OJ 

1 

u 

1 

OJ 



{ff 



D 



-1 



[B] [A,]-'{f} 



[B] [A,]-'{f} 



- [Ad] + [B] [A,]-^ [B] 



1 


*- 







{/r[A,]-MB] -[Ad] + ^[B][AJ-Mb] 



{/} 



[Ad] 



-^[B] [A,]-^ [B] 



\B]-' [A. 



-1 



{/} 



Ad] iB]-VAJ + ^[B 



-1 



{/} 



(36) 



and, 



^2 



1 
1 
1 



iff 

{fV 

iff 



[A,] + ^ [B] [Ad]-Mb]] '[B][Ad]-U/} 



[Ad] [B]-^ 



[A,] + - [B] [Ad 



-1 



B] 



-1 



{/} 



n -1 



[Ad] [B]-MA,] + -[B] 



{/} 



(37) 



which show that SI = S2. It can be seen that the matrices [Ad], [Ap], and 
[B] are Hermitian regardless of the asymmetric nature of the unit celL While 
asymmetries related to dimensions affect the off diagonal terms in a way which 
preserves the Hermitian nature of the matrices, asymmetries related to ma- 
terial properties only affect the diagonal terms, thereby, preserving the Her- 
mitian characteristics. This in effect means that, for any asymmetrical unit 
cell, D and p are always real-valued and that the coupling terms in (13T| [32!) 
are always complex conjugates of each other. This also ensures that the total 
energy given by, 
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E = -[{ar{e) + {a){er + {ur{p) + {u){pr] 



(3J 



2D{a){a)* + 2p{u){u)* + 2S2{(t){u)* + 25i(a)*(it) 
2D (a) {ay + 2p{u){uy + 2S*((t)(m)* + 2Si{ay{u) 



D{a){ay + p{u){uy+ S,{ay{u) +S^{ay{u) 



is always real. Now, since E is also non-negative, and that the average stress 
and velocity can be prescribed essentially arbitrarily and independently, it 
follows that D and p are also positive. 



2.2 Bloch-form Overall Spatial Variation and Dispersion Relations 



We now consider a special case of an infinite homogenized elastic solid with a 
layered microstructure, and seek conditions under which it supports periodic 
waves of the Bloch-form spatial variation of the following form: 



(a)(x) = (a)e^^- {u)ix) = {u)e^'>- ^^^^ 
{e){x) = {e)e'^^; {p){x) = (p)e'^^ 



The overall field equations then become. 



which yield, 



q{a) = -uj{p)- q{u) = -u{e) (41) 

From these equations and the overall constitutive relations fl?I| 15^ it now 
follows that, for this homogenized continuum to admit solutions of the form 
( 139|) . the following dispersion relation must necessarily be satisfied: 



f(A ^ (1 + VpSi){l + VpS2) .^2) 

\q) P Dp 

where Vp is the phase velocity. The above equation is used to evaluate the 
dispersion relation for a layered composite by substituting from Eq. (1331) . 
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The real frequency-wavenumber pairs are then used to evaluate the effective 
parameters of Eq. ( 13T| [32]) and these parameters would have the properties 
discussed in the preceding subsection. 

If we wish, as in [5], to directly relate the overall strain, (e), to the overall stress, 
(cr), and the overall momentum, (p), to the overall velocity, {ii), according to 

(e) 
(P) 

then we have. 



^eff 



and they do satisfy the same dispersion relation Eq. ( I42l) . now rewritten as. 

In this representation however, the effective parameters defined by fl43p are 
real only for symmetric (elastic) unit cells and become complex-valued for 
asymmetric unit cells. Also the parameters D'^^ and p'^^ are obtained using 
the special form of the field equations f HTj) while the constitutive relations fl3T| 
15^ have broader applicability. 

A significant advantage of our homogenization technique is that the point- 
wise solution of the field equations is not required a priori, and that, this 
solution can be obtained as a by product of the approach. In fact, within each 
subregion, the value of each field variable is given by adding to each periodic 
part in Eqs. fl2 ^[21|) . the corresponding average value. We do not need explicit 
expressions for the mode-shapes, but rather can extract these from the final 
results, if desired. The dispersion relations are obtained by substituting Eqs. 
( 13T| |32|) into Eq. (H2|) . These relations then relate the frequency, u, and the 
wavenumber, q, for the special Bloch-form of the periodic waves, as exemplified 
in what follows. 



D 



1 + VpSi 



P 



1 + Vr,S2 



(43) 



(44) 



(45) 
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2.3 Example: Dispersion calculation for a layered composite 



To illustrate dispersion calculations by our micro-mechanical formulation, con- 
sider the case of a layered composite (Fig. [T]) with harmonic waves traveling 
perpendicular to the layers. 



Direction of 
wave propagation 



X 



Material 1: p^lOOOkg/m' 
E=2 X 109 Pa 
Thickness=5mm 

Material 2: p=3000 kg/m' 
E=200 X 109 Pa 
Thickness=5mm 



Fig. 1. Schematic of a bi-layered composite 

The exact dispersion relation for longitudinal waves in a layered composite 
was given by Rytov [5]. Here we compare the dispersion curves calculated 
using our micro-mechanical formulation with those of the exact solution. A 
unit cell, for this case, consists of 1 layer of Material 1 and 1 layer of Material 
2. We divide the unit cell into N = 30 subregions (15 divisions per layer) and 
use constant values for eigenstress and eigenvelocity within each subregion. 

2 2.00E+06 

>: 

^ 1.50E+06 

(U 
3 

g" l.OOE+06 

iS 5.00E+05 
3 
u 

i— 

U O.OOE+00 

0.5 1 1.5 2 2.5 3 

Normalized wavenumber 

A Micromechanical Formulation # Exact Solution 

Fig. 2. Dispersion curve comparison of the Micro-mechanical formulation with the 
exact results 

The first three modes in the dispersion curve are compared in Fig. ([2]). It can 
be seen that the micro-mechanical formulation gives very acceptable results 
for the first three modes when appropriate discretization is used. There are 
indications that the results of the micro- mechanical formulation start diverging 



«■■■ 1 



-m • • 
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from the exact results as one considers higher frequencies. This is evidenced by 
the increasing difference in the results as one moves up the third branch. This 
is of course an issue of discretization and the micro-mechanical formulation 
is expected to give accurate dispersion results at any frequency given a high 
enough discretization. 

It should be mentioned here that exact dispersion relations are only available 
for fairly simple geometries like layered composites. The lack of exact solutions 
for more complex cases, therefore, necessitates accurate and efficient numerical 
schemes. In a recent paper [5], a mixed method is used to obtain the point- wise 
values of the field variables and to formulate the corresponding integration- 
based homogenization of layered composites. This method was shown to be in 
good agreement with homogenization based upon the integration of field vari- 
ables as calculated from the exact solution. Although homogenization based 
on the mixed variational method provides an accurate and efficient numer- 
ical alternative to using the exact solution, the field variables still need to 
be evaluated at multiple points within the unit cell. Here we show that the 
micro-mechanical approach to homogenization, in addition to precluding the 
need for point-wise evaluation of the field variables, also converges to homog- 
enization results based on both the mixed variational method and the exact 
solution; a brief summary of integration of field variable based homogenization 
is presented in Appendix |X1 



3 Constitutive relations of layered composites 



To illustrate the micro-structural method of calculating the constitutive pa- 
rameters of homogenized composites, examples of 4-layered symmetric and 
asymmetric composites are considered. These examples also serve to show 
that homogenization based on the micro-structural formulation converges to 
the results of homogenization based on the integration of field variables using 
the exact solution and/or using the results of the mixed variational method. 



3.1 -^-layered symmetric composite: Comparison with homogenization based 
on the mixed variational formulation 



We consider a 4-layered composite in order to show that the micro-structural 
homogenization gives results which converge to homogenization based on the 
integration of field variables as derived from the mixed variational formula- 
tion. We stressed that for this integration-based homogenization, the effective 
parameters are calculated from Eqs. ( lA.lOl IA.11|) using the mode-shapes cal- 
culated from the approximate mixed variational formulation of [TB], whereas 
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for the micro-structural homogenization the constitutive parameters are cal- 
culated from Eqs. fH^ |45!) directly. It should be noted that while the effective 
parameters calculated from Eqs. (lA.lOl lA.lip require integration of the field 
variables through the unit cell, effective parameters calculated from Eqs. fl44[ 
H5|) do not require any integration. For symmetric unit cells such as Fig. ([3^), 
Eqs. dm |45|) result in real valued effective parameters and produce the same 
constitutive relations flA.10tlX.ll|) reported in literature [5]. Expressions flA.lOt 
lA.lip however do not clearly display the coupling between average strain and 
average velocity, and average momentum and average stress, which is non-zero 
even in the symmetric case ([231) • The micro-structural approach directly pro- 
duces the effective parameters of Eq. (!33l) without any a priori assumption 
about the structure of the constitutive relations. These effective parameters 
can always be used to calculate and p'^^, if desired. 



IIJ 



Material 1: Density=1000 kg/nr^ 

Compliance=1/(8X10'' )Pa''' 



Material 2: Density=300 kg/nr' 

Compliance=1/(3X1 0^)Pa"'' 

Material 3; Density=8000 kg/m^ 

Compliance=1/(300X10'''')Pa"'' 




0.5 1 1.5 2 2.5 3 

Normalized Wavenumber (Q) 



Fig. 3. 4-layered symmetric composite: a. Schematic of a unit cell; b. Dispersion 
curve for the first 2 propagating branches 



Fig. ([3^) shows one unit cell of the 4-layered composite under consideration. 
It is composed of a heavy and stiff layer sandwiched between two layers of a 
light and compliant material with the whole assembly being embedded in a 
heavy and stiff matrix. Here we calculate the effective dynamic properties of 
the composite using the micro-structural method dlU H5|l with increasingly 
more refined discretization and compare the corresponding homogenization 
results with those obtained from Eqs. flA.10[ lA.lip . 



Fig. dl^) shows the homogenization results for the first two propagating branches 
calculated from the integration of field variables as derived from the mixed 
variational formulation. Fig. (IDd) shows the homogenization results from the 
micro-structural formulation. This example also serves to show that the ini- 
tial choice of the homogenizing parameters (po, -Do) is immaterial to the fi- 
nal homogenization results. In this example po and Dq are taken to be the 
average properties of layers 1 and 2. Therefore, the discretization of hetero- 
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Fig. 4. Effective Parameters comparison: a. Homogenization results from field in- 
tegral method as applied to the mixed variational formulation; b. Homogenization 
results from the micro-structural method with increasingly more refined discretiza- 
tion 

geneity is applied to the entire unit cell and every layer is discretized. The 
effect of increasingly more refined discretization is shown in Fig. (jib)- The 
number of terms in the Fourier expansion (^) is kept constant at 10 for each 
level of discretization. From Fig. (jib), it can be seen that when each layer 
is considered as just one subregion of constant eigenstress and eigenvelocity 
(A^ = [1, 1, 1, 1, 1]), the homogenization results deviate significantly from the 
results of the field integration method. As finer discretization is introduced in 
each layer (A^ = [3,3,3,3,3]), homogenization based on the micro-structural 
method tends to converge to the results of the integration method using the 
mixed variational formulation. For a suitably fine discretization, especially for 
the case when the more compliant layer is modeled with greater number of 
subregions (A^ = [5,15,5,15,5]), the results of the micro-structural method 
approach those of the mixed formulation, with any desired accuracy. 



3.2 J^-layered composite with an asymmetric unit cell 



We now present an asymmetric unit cell which serves to illustrate the nature 
of the coupling parameters in Eq. (133|) . For a symmetric unit cell, all effective 
parameters calculated from Eq. ( I33l) are real-valued and yield real-valued D^^ 
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and p^^. It is for the asymmetric unit cell that the coupling terms in Eq. ( 133|) 
become complex- valued resulting in complex- valued D^^ and p'^^. 

Fig. (E^) shows one unit cell of the 4-layered asymmetric composite under con- 
sideration. It is composed of a heavy and stiff layer sandwiched between two 
equally thick layers made of different compliant materials, with this assembly 
being embedded in a heavy and stiff matrix. Here we calculate the effective dy- 
namic properties of the composite using the micro-structural method (1441145 p . 



a b 




Normalized Wavenumber (Q) 



Fig. 5. 4-layered asymmetric composite: a. Schematic of a unit cell; b. Dispersion 
curve for the first 2 propagating branches 

Fig. (E]) shows the homogenization results for the first two propagating branched 
calculated from the micro-structural formulation (Eq. |33|) . In this example po 
and Dq are taken to be the properties of layer 1. Therefore, the discretization 
of heterogeneity applies to the central three layers. The number of terms in 
the Fourier expansion (^) is kept constant at 10 and the discretization level 
of the central three layers is (A^ = [15, 10, 15]). The real and imaginary parts 
of the calculated parameters are indicated with different symbols. Figs. (l6K,d) 
clearly show that D and p are purely real for all propagating frequencies and 
Figs. (I6b,c) show that Si and 5^2 are complex conjugates of each other. 



4 Conclusions 

A homogenization method based on micromechanical considerations [T|TB] is 
presented. This method does not require the point-wise values of the elasto- 
dynamic field variables within a unit cell, instead, it produces those values 
to any desired degree of accuracy, if needed. Using the local conservation and 
kinematical relations together with the local constitutive relations, we have 
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Fig. 6. Effective parameters calculated from the micro-structural formulation, a. Z?, 
b. 5i, c. 5*2, d. p. 



systematically deduced the overall constitutive relations for the homogenized 
elastic solid that depend on the frequency and wavenumber as two indepen- 
dent parameters. These parameters can be related to one another in solving 
any special overall initial-boundary-value problems associated with the ho- 
mogenized continuum. For example, for a Bloch-form wave propagating in an 
infinite homogenized continuum (microstructurally layered periodic compos- 
ite), one arrives at the dispersion relations that connect the wavenumber to 
the wave frequency, as necessary conditions for the existence of such solutions. 
In general, for elastic waves in a layered unit cell, the average strain and aver- 
age momentum are given in terms of the average stress and average velocity 
by a pair of coupled constitutive relations. We have proved that the overall 
compliance and mass-density in these relations are always real-valued, and 
that the corresponding coupling terms are always complex conjugates of each 
other. For a symmetric unit cell, the coupling terms are real- valued and equal 
to one another. 

For illustration, we have used the overall constitutive relations to calculate 
the dispersion relation of a layered composite and have shown that the results 
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converge to the exact solution for the case of a 2-layered composite. For the 
special problem of a Bloch-form wave propagating in an infinite homogenized 
solid, one may use the overall conservation and kinematical relations to de- 
fine another set of effective parameters, D*^^ and p'^^, that directly relate the 
average strain to the average stress, and the average momentum to the aver- 
age velocity, similar to the corresponding quasi-static case; see Nemat-Nasser 
and Hori [1] . We have shown that these parameters are complex- valued for any 
non-symmetric unit cell. We have presented a symmetric 4-layered example 
showing that these parameters as calculated from the micromechanical method 
converge to those calculated by the integration of field variable approach [5] . 
The microstructural method we presented here however, may be confidently 
used to homogenize unit cell of more complex microstructures in 2- and 3- 
dimensional composites where the exact dispersion relations are not known 
and the point-wise solution to the field equations is not available, as we shall 
show in a subsequent work. Finally, we presented an asymmetric 4-layered ex- 
ample where we have shown that the coupling terms are complex conjugates 
of each other and D and p are real-valued, whereas the corresponding D^^ 
and p'^^ are complex-valued. 
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A Homogenization by integration of field variables 

Nemat-Nasser et. al. [5] proposed a homogenization method of periodic elastic 
composites based upon the integration of field variables. Here we present a 
brief summary of the basic equation. For harmonic waves traveling in a layered 
composite with a periodic unit cell Q = {x : —a/2 < x < a/2} the field 
variables (displacement, velocity, stress, and momentum) take the Bloch form 
given by Eq. The dynamic equilibrium gives, 

Vcr + iup = (A.l) 
where V denotes differentiation with respect to x. 

A.l Effective properties 

Multiply Eq. (lATj) by e''"^ and use Eq. to obtain, 

V (a(x)e^«(^-^)) + iwp(x)e^«(^-^) = (A.2) 

Introduce the change of variable y = x — X and average with respect to X 
over a unit cell to obtain, 

((f^)e"''') + ito{p)e"''^ = (A.3) 

where 

1 r+a/2 I r+a/2 

(a) = - a{x)dx; (p) = - p{x)dx (A.4) 

a J -a/2 a J -a/2 
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Now define the mean stress and mean momentum density as, 



(a)(a;) = (a)e^''^; = (p)e^^^ (A.5) 

Observe that the mean stress and mean momentum density satisfy exactly the 
overall equation of motion, 



V{a) + iuj{p) = (A.6) 
Define the effective mean displacement as. 



= (M)e''?"; {u) = - r"^\{x)dx (A.7) 



a J-a/2 

Then the effective mean strain and velocity are given by, 

{i){x) = iq{u){x); {u){x) = —iu{u){x) (A. 8) 

Finally, the averaged stress-strain relation becomes, 

D^^{a){x) = {e){x) (A.9) 
Based on Eq. flA.5l lA.Sp . D"^ is given by, 

D^^ = ^^ (A. 10) 

Similarly, the effective density (p*^^) is given by, 

—IU!{U} 

Since (a) and (p) satisfy the equation of motion (Eq. lA.Gp . it is seen that the 
effective material parameters defined above automatically satisfy the disper- 
sion relation, 
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